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We present the performance of searches for gravitational wave bursts associated with external 
astrophysical triggers as a function of the search sky region. We discuss both the case of Gaussian 
noise and real noise of gravitational wave detectors for arbitrary detector networks. We demonstrate 
the ability to reach Gaussian limited sensitivity in real non-Gaussian data, and show the conditions 
required to attain it. We find that a single sky position search is ~20% more sensitive than an 
all-sky search of the same data. 
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I. INTRODUCTION 

Searches for transient gravitational waves (bursts) typ- 
ically fall into one of two categories: all-sky untriggered 
searches, which scan the entire sky and search through- 
out the available data; and triggered or directed searches, 
which historically analyze only a single point on the sky 
corresponding to some astrophysical source of interest 
(see for example Refs. [3, 0]). However, some instru- 
ments provide only an approximate sky location for ex- 
ternal triggers, which requires scanning a relatively large 
patch on the sky for an associated gravitational wave. For 
instance, this is the case for gamma-ray bursts (GRBs) 
localized by the Gamma-ray Burst Monitor (GBM) [g, |4| 
on Fermi, and for some high energy neutrino candi- 
dates 0,3. 

We present an implementation of a gravitational wave 
burst (GWB) search which is able to scan arbitrary sky 
patches, and demonstrate how the sensitivity of the anal- 
ysis varies with the sky region searched. In Gaussian 
noise we find that the sensitivity is a function of the 
total signal-to-noise-ratio received by the gravitational 
wave detector network, which is the expected result for 
an optimal search in Gaussian noise. However for real 
non-Gaussian gravitational wave detector noise, we find 
a different dependence of the sensitivity, due to the re- 
quirement that the gravitational wave signal be seen by 
at least two detectors to be distinguished from spurious 
noise transients. We obtain an empirical formula for the 
sensitivity of the search in real noise as a function of 
the sky position for arbitrary detector networks. In par- 
ticular, a sensitivity loss of ~20% is observed between 
searching at a single sky position and an all-sky search 
of the same data. 



Finally, we demonstrate the ability to reach Gaussian 
limited sensitivity in real non-Gaussian data, and show 
the conditions required to attain it. Specifically, we find 
that all of the detectors present in the network need to 
have comparable sensitivity. Moreover, we show that in 
practice adding a detector to a gravitational wave detec- 
tor network may actually reduce the search sensitivity for 
some sky areas when analyzing real data. 

We begin in Sec. |TT] with a brief introduction to co- 
herent searches for GWBs. We follow in Sec. Mil with 
details of the detection statistic and sky scanning algo- 
rithm used in this paper. The spurious-noise rejection 
tests are presented in Sec. IIVI In Sec. [V] we present the 
performance of the search for GWBs; in particular, the 
sky dependence is discussed in Sec. IV A) the comparison 
between real and Gaussian noise is shown in Sec. IV Bl 
and the dependence on the size of the sky region is pre- 
sented in Sec. IV CI We conclude with some comments on 
the implications of these results in Sec. IVII 



II. COHERENT ANALYSIS OVERVIEW 

Coherent analysis of gravitational wave data was orig- 
inally introduced in Ref. Q- Since then it has proven 
to be an effective method for GWB searches in LIGO- 
Virgo data, and it is now the dominant search method- 
ology [l|, [H, 0, l8rfl3j . Here we give a short overview of 
coherent analysis in order to introduce notation used in 
the following sections. 

For a gravitational wave (h + (t), h x (t)) incoming from 
a sky location Cl the calibrated data from a gravitational 
wave detector a are of the form 
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d a (t + At a (fi)) = 
F+(Cl)h+{t) + F* (£l)h x (t) + n a (t + At a (Cl)) . (1) 
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Here (f2) , F* (Tl) are the antenna response functions 
of the given detector to the plus (+) and cross (x) polar- 
ized gravitational waves, n a is a time series of detector 
noise, and At a (£l) is the gravitational wave travel time 
between the position r a of the detector and an arbitrary 
reference point rb: 



At a (£l) = ~(r -r a )-£l. 
c 



(2) 



For the case where the sky location f2 is known a priori, 
the first step of the analysis is to time shift the data d a 
by the known Ai Q (f2) in order to obtain a gravitational 
wave contribution which is synchronous between the dif- 
ferent data time series. The data from each detector 
are then whitened and decomposed in a time-frequency 
representation, e.g. using a short Fourier transform or 
a wavelet decomposition, where the time-frequency ba- 
sis functions typically have length between several mil- 
liseconds and several hundred milliseconds. For a given 
time-frequency pixel (basis function) of center time t and 
center frequency / the obtained decomposition can be 
written compactly as 



d = F + /l4 



F x h> 



n . 



(3) 



where boldface symbols denote the vector of whitened, 
time-frequency decomposed time series in the D dimen- 
sional space of detectors: 
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Here Si(f) is the one-sided noise power spectrum of de- 
tector i. Note that the gravitational wave contributions 
h +1 h x are the projection on a time- frequency basis func- 
tion without any whitening. 

The basis used in Eqs. (|4]) to describe vectors in the 
D dimensional space of detectors is not adapted to the 
gravitational wave contribution. Given that the data are 
whitened and detector noise can be assumed to be uncor- 
rected between detectors, the vector n has an identity 
covariance matrix, which is invariant under change of or- 
thonormal basis. Hence we can construct a new adapted 
orthonormal basis, in which the first two vectors span the 
gravitational wave plane [13] generated by the F + and 
F x vectors, and the D — 2 remaining vectors span the 
null space, the space orthogonal to the gravitational wave 
plane. 



This basis can be further refined in various ways, such 
as by choosing the first two vectors along the directions 
of maximal and minimal response to a linearly polarized 
gravitational wave. These two directions are orthogonal 



and correspond to the dominant polarization choice 15 1 
of the arbitrary gravitational wave polarization angle ref- 
erence, which is a particular choice of the definition of the 
plus and cross polarizations. We denote the antenna re- 
sponse vectors for this special polarization choice by / + 
and / x . They have the properties 



i/ + i 2 >i/ x r 

/+•/* =o. 



(5a) 
(5b) 



Note that the choice of which of the plus and cross vector 
has a larger amplitude is purely conventional. The unit 
vectors of our adapted basis are hence e + = f + /\f + \ 
and e x = / x /|/ x |, complemented by vectors spanning 
the null space, for instance e" = e + A e x for the case of 
3 non-aligned detectors. 

An alternative basis choice may be appropriate when 
we have prior information on the expected gravitational 
wave polarization. For example, for a circularly polar- 
ized gravitational wave signal the projection onto a time- 
frequency basis function of the two polarizations are re- 
lated by 



±ih 



+ ! 



(6) 



depending on whether the signal is left or right-hand po- 
larized. Hence in the detector space the gravitational 
wave contribution will lie along either the left or right- 
handed response vectors 



with corresponding unit vectors e 
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(7) 
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To construct a left or right-handed basis we 
use the vectors orthogonal to the response vectors in the 
gravitational wave plane 



r 



nO 



\r 



. (8) 
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with the corresponding unit vectors e n(J 
and e"° = f n0 /\f n0 \. 

Either basis (dominant or circular) can be used to con- 
struct two types of statistics: detection statistics and co- 
herent consistency statistics. A detection statistic is used 
to ranks events as more consistent with a given model 
of signal buried in noise than a model of noise only. It 
is usually expressed as the log-likelihood ratio between 
these two models. Coherent consistency tests, on the 
other hand, are used to reject spurious noise transients 
which are usually not included in the noise model of the 
detection statistic. We discuss the formulation of detec- 
tion statistics in Sec. IIIIl and coherent consistency tests 
in Sec. IIV1 for the moment we note that these statistics 
are only functions of the data vector d and the antenna 
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response vectors / + and / x or f^,f ( ^ and Z™ ,/"^ 1 . 
The value of each statistic is computed independently 
for each time-frequency pixel, and the resulting values of 
the detection statistic over the array of time-frequency 
pixels is used to define gravitational wave events. Specif- 
ically, all time-frequency pixels with detection statistic 
above a certain threshold are clustered to form events, 
for instance using nearest-neighbor clustering 16]. The 
final detection statistic of such an event is simply the sum 
of the detection statistic of all the pixels composing the 
event due to the additive properties of the log-likelihood 
ratio. The other statistics are also summed over the clus- 
ter of pixels composing the event. 

This event generation procedure is used on numerous 
background samples (generated from real data using the 
time slide technique) and signal samples (generated by 
adding simulated gravitational wave signals to real data). 
The obtained background and signal events are used to 
tune the coherent consistency tests to reject the tail of 
spurious noise transients inconsistent with the gravita- 
tional wave signal hypothesis. Independent samples of 
background events are then used to estimate the distri- 
bution of background events that survive the consistency 
tests, which is used in turn to define the statistical signifi- 
cance of any candidate gravitational wave events from the 
analysis of data coincident with the external astrophys- 
ical trigger. An independent sample of simulated signal 
events is used to estimate the sensitivity of the analysis 
as a function of gravitational wave signal amplitude, and 
to construct upper limits on the gravitational wave sig- 
nal amplitude whenever no significant gravitational wave 
event is found. 



III. DETECTION STATISTIC 

In a gravitational wave search, the detection statistic 
is used to ranks events as more consistent with a given 
model of signal buried in noise than with a model of noise 
only. The detection statistic is often based on some mea- 
sure of the energy in the data, motivated by a likelihood- 
ratio analysis. For the present analysis, we construct 
a detection statistic following the Bayesian formalism 
of Ref. [l7j . The signal model is a circularly polarized 
gravitational wave signal with Gaussian amplitude dis- 
tribution of width er/j, and the noise model is Gaussian. 
The circular polarization assumption is well motivated 
for some astrophysical sources, for instance gamma-ray 
bursts, as discussed in Sec. IIV Bl For a right circularly 
polarized signal we obtain the log-likelihood ratio 

2L(d\ 0,a h ) = I 6 "' d [ 2 -log(l + <7g|/°| 2 ) , (9) 

and the left circular polarization log-likelihood ratio 
L(d\ 0,(7/0 has an analogous form. The final likelihood 
ratio is obtained by marginalizing over the left versus 
right choice, and over a discrete set A of ah covering the 



range [10 23 ,10 21 ] Hz 1 / 2 of realistic detectable signals. 
The detection statistic used thus has the form 

•^detection = L[d) = 

log 2\A\ i ex P L ( d \ 0,(T h ) + expL(d\ <3,a h )}. (10) 

cr h eA 

The discussion so far has assumed we know the sky 
position of the gravitational wave source a priori. For 
some searches this is indeed the case, such as for gamma- 
ray bursts detected by the Swift satellite [Hj]. In other 
cases, such as untriggered all-sky searches 0, l^-Hlj. Cl 
is not known, or may only be constrained to some large 
region of the sky. An example of the latter is gamma-ray 
bursts detected by the GBM, which has relatively large 
sky location systematic uncertainties of a few degrees [3[ 
and statistical errors of up to ~ 10 degrees depending on 
the 7-ray flux and spectrum. An error of this size in the 
sky location used to synchronize the data time series from 
gravitational wave detectors causes timing discrepancies 
of up to several milliseconds. A potential gravitational 
wave signal at a few hundred Hertz could easily be shifted 
by a quarter of a period or more between a pair of gravi- 
tational wave detectors, and the signal could be rejected 
by a coherent consistency test. 

The standard solution in gravitational wave coherent 
searches is to repeat the analysis over a discrete grid of 
sky positions covering most of the source sky location 
probability distribution. Here we use a simple regular 
grid, composed of concentric circles around the best es- 
timate of the source sky location, which covers at least 
95% of the sky location probability distribution. The 
constant grid step is chosen so that the timing synchro- 
nization error between any sky location in the error box 
and the nearest analysis grid point is less than 10% of 
the period for the highest frequency gravitational wave 
signals included in the search. This is small enough to 
limit the amplitude SNR loss due to timing error to be 
less than 10%. 

Gravitational wave triggers are produced indepen- 
dently for each sky position grid point. The detection 
statistic ^detection, written as a log-likelihood ratio be- 
tween a signal and a background model, is penalized by 
the probability pem(^) of the trial sky location being the 
true one, by adding the logarithm of that probability to 
the detection statistic 

^penalized = ^detection + logPEM^J ■ (H) 

The reconstructed source sky position for a given signal 
is the sky position for which the trigger has the largest 
penalized detection statistic. Only that maximal trigger 
is kept by the analysis of the grid of sky positions. 

As a simple model of 7-ray satellite errors we use a 
Fisher probability [l9[ 

PEM(fi) = PFi S her(0; k) = ~ e * C ° S 9 , (12) 

£jrv <2l — n. 
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where 6 is the angle between the best estimate sky loca- 
tion and the analyzed one. The parameter n is chosen so 
that the 95% coverage radius of this Fisher distribution 
is equal to the 95% coverage radius for a given GRB sky 
location reconstruction (with statistical and systematic 
errors added in quadrature). This model is a reasonable 
approximation for localization performed by a single 7- 
ray spacecraft, such as Fermi or Swift, however it may not 
apply to other instruments, for instance to localization by 
the Third Interplanetary Network of satellites [20| • 

IV. COHERENT CONSISTENCY TESTS 
A. General framework 



extended to projection vectors that are not in the null 
space, but in the gravitational wave plane. The incoher- 
ent energy will remain comparable to the coherent energy 
for the case of noise transients, but for gravitational wave 
signals it will be much smaller or larger than the coherent 
energy. The previously proposed extension [25| is to use 
the plus and cross polarization directions in the dominant 
polarization frame. Along e + a coherent buildup of en- 
ergy is expected for most gravitational wave signals; this 
projection corresponds to the hard constraint introduced 
in Ref. [jq . On the other hand, for many network config- 
urations and large fractions of the sky |/ x | <C |/ + | [15j . 
and the projection on the e x vector can be considered as 
an effective null stream. 



Detection statistics such as Eq. (|10[) are usually con- 
structed assuming the background detector noise is Gaus- 
sian. These statistics do not take into account spurious 
noise transients, mainly because no good model of these 
transients is available. However these transients are un- 
correlated between the different gravitational wave de- 
tectors, and powerful coherent consistency tests to reject 
them can be constructed on that basis [2ll - l24j |. 

An effective method for rejecting noise transients is 
to project the gravitational wave data vector d onto the 
null space, and to compare the squared magnitude of this 
projection with the autocorrelation terms of that projec- 
tion [23] . For simplicity let us consider the case of a one 
dimensional null space along a vector e". The squared 
magnitude of the projection on this vector, also called 
the coherent null energy, is 

£ n =|e"-d| 2 = £Ce£OV (13) 

The autocorrelation part of the null energy, called the 
incoherent null energy, is 

4 = E«^a = El e S| 2 K| 2 ' ( 14 ) 
a a 

For a strong gravitational wave signal the contribution 
from the different detectors cancel each other in the co- 
herent null energy, but stay present in the incoherent null 
energy, hence E n <C I n is expected. For noise transients 
the contributions of each detector are expected to be un- 
correlated, hence E n ~ /„. Thus, a threshold on the 
ratio of the incoherent and coherent null energy can be 
used to reject noise transients. 

Several extensions of this framework have been previ- 
ously implemented and discussed in Ref. [25j . First, the 
separation between the coherent and incoherent energy 
for both noise transients and gravitational wave signals 
depends on the value of the incoherent energy, that is 
the strength of the deviation from the Gaussian noise 
hypothesis, hence a more complicated separation line in 
the coherent/incoherent energy plane than a simple ratio 
is used in practice. Second, the framework has also been 



B. Extension to circular polarization 

The main issue with the previously described projec- 
tions is that none of them is effective for the case of two 
gravitational wave detectors which see roughly indepen- 
dent linear polarizations [26], which occurs frequently for 
networks consisting of one LIGO detector plus Virgo. In 
that case the coherent consistency tests based on the 
plus and cross energies perform poorly, as the cross- 
correlation terms are small, and the incoherent and co- 
herent energy are roughly equal for both gravitational 
waves and noise transients. This issue was noted in the 
all-sky GWB search of 2005-2007 LIGO- Virgo data [J, 
and can also be seen in the poor upper limits for the 
search in association with GRBs of the same data 0] ■ 

Here we propose to expand the possible projection vec- 
tors using the assumption of a circularly polarized grav- 
itational wave signal. This imposes correlations between 
otherwise independent linear polarizations, and allows us 
to construct effective consistency tests even for the case 
of two strongly misaligned gravitational wave detectors. 
Also, we note that the circular polarization assumption is 
well motivated for many astrophysical scenarios. For in- 
stance, for gamma-ray bursts the expected beamed emis- 
sion of the gamma rays means that the progenitor is seen 
roughly along its axis of rotation. Furthermore, gravita- 
tional wave emission models which could be seen at extra- 
galactic distances emit predominantly circularly polar- 
ized gravitational waves along their rotation axis [27H3l| . 

We consider the projection onto the manifold of circu- 
larly polarized gravitational waves, which is formed by 
the two complex lines along the left and right handed 
polarization directions. The magnitude of the projection 
is simply the maximum of the projection on the left and 
right handed response unit vectors and , which we 
defined in Sec. [Hi 

Circular = max (|e° • d\ 2 , \e° ■ d\ 2 ) . (15) 

This projection can be compared to its incoherent part 

^circular = = i (16) 
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FIG. 1. Strain noise spectra from the LIGO detectors used 
in this study. Note that H2 data are treated as if the de- 
tector was located at the Virgo site in order to simulate the 
gravitational wave detector network operated in 2009-2010. 

which is equal for the right and left projections. Using 
the pair of incoherent / coherent energies ^circular 5 ^circular 
a consistency test can be constructed analogously to the 
ones based on the plus or cross energies. 

Furthermore, given the circular polarization assump- 
tion a null test can be constructed by considering the 
unit vectors e"° and e nQ that are orthogonal to the cir- 
cular projection and defined in Sec. [TTJ The null circular 
energy is defined as the minimum of the magnitude of 
the projection onto these two unit vectors 

Kull circular = min (|e"° ■ d\\ |e™° ■ d\ 2 ) . (17) 

As in the previous cases an incoherent counterpart can 
be defined, with the autocorrelation terms of the null left 
and right projection being equal. A consistency test can 
then be constructed using E nua c i, cu i ar , I uu n circular- 

V. ANALYSIS PERFORMANCE 

The analysis methodology described above has been 
implemented in X-Pipeline 25;], a software package de- 
signed for GWB searches in association with external 
astrophysical triggers. Earlier versions of this package 
have been used to search for GWBs associated with 
GRBs @, H] . 

In order to understand the sky dependence of the sen- 
sitivity of GWB searches, we study the performance of 
X-Pipeline in searching for circularly polarized gravi- 
tational wave bursts. Specifically, we characterize the 
performance for a typical GRB-trigger scenario for the 
most recent science run of the LIGO- Virgo network (32r - 
l33 | . from 2009-2010. Since data from that period have 
not yet been released for data analysis performance stud- 
ies of the type presented here, we use as a proxy a 3 hour 



long sample of LIGO data from 23 February 2006. At 
that time the two detectors at the Hanford site and the 
detector at the Livingston site were taking science qual- 
ity data; the spectra of the data at the center of that 
sample are shown on Fig. [TJ To be representative of the 
full network of large scale interferometric gravitational 
wave detectors which was operational during 2009-2010, 
we use the data from 2 km detector Hanford site as if that 
detector was located at the Virgo site, hence throughout 
this article we will denote these data as from VI. The 
factor ^2 difference in sensitivity is roughly representa- 
tive in the difference in sensitivity between Virgo and the 
two 4 km LIGO detectors (HI and LI) in the 2009-2010 
data set. 

These data are used to generate background data sam- 
ples using the time slide method, and to generate simu- 
lated signal samples by adding circularly polarized sine- 
Gaussian and compact binary inspiral waveforms into 
these data. We perform full autonomous analyses as used 
in real externally triggered gravitational wave searches 
to determine the sensitivity for these signal models. We 
use a 660 s long time window around fiducial external 
triggers, and search for gravitational waves over the fre- 
quency band 64—500 Hz. These are the parameters 
used in the search for gravitational waves associated with 
GRBs in 2009-2010 LIGO- Virgo data QJ]. For reference, 
we also use simulated Gaussian noise with the same spec- 
tral properties. With this reference the effect of non- 
Gaussian noise transients present in real data can be as- 
sessed. 



A. Single sky position analysis 

To study the effect of different antenna pattern con- 
figurations, we perform the analysis of fiducial external 
triggers well localized to different points in the sky, which 
is representative of gamma-ray bursts reported by Swift. 
We use the same time for all these external triggers, but 
select sky locations which probe different combinations 
of contributions from the various detectors in the grav- 
itational wave detector network. The fiducial sky loca- 
tions and gravitational wave detector networks used in 
this study are shown on the 3 panels of Fig. O 

In order to have an a priori measure of the signal 
strength with regard to the detectors noise, we define 
the coherent network SNR 35] 

p° = J E pi ( 18 ) 

y ct E network 

of a gravitational wave signal model as the root-sum- 
square of the individual SNRs 

/=4 r i^M/)+j?M/)i 2 d/ (19) 

Pa ~ Jo s a (f) 1 1 ] 

of that signal in each detector a across the network. This 
figure of merit can be used as a predictor of the gravi- 
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FIG. 2. Sky dependence of the penalty factor A c /p c [Eq. (|22|l ] 
assuming that VI is a factor 2 less sensitive than HI and LI. 
From top to bottom the penalty factor for respectively the 
H1L1, L1V1 and H1L1V1 networks are shown. The coordi- 
nates used are the longitude and latitude in a Mollweid pro- 
jection, as the detectors are fixed to the Earth. The location 
of the detectors in each network are marked by blue dots; 
for each detector the projection on the celestial sphere of the 
location (and its opposite) are the points of maximum an- 
tenna response. For the 2 detector networks the zero points 
of the penalty factor correspond to the null of the antenna 
response for one of the two detectors. The plus marks show 
the sky location of the fiducial triggers analyzed to produce 
Fig. [3] whereas the circles show the large sky regions used for 
the study shown in Fig. [S] We omit the map for the H1V1 
network as it is very similar to the LI VI map. 



tational wave signal amplitude (distance to the source) 
needed for the signal to be detected. 

We define the detection sensitivity of the search for a 
particular signal model as the distance d 50 o/ at which 
that signal is found with 50% efficiency while holding 
the background false-alarm probability fixed at 1% per 
on-source window (in this case 1% per 660 s, or a false 
alarm rate of 1.5 x 10~ 5 Hz). This detection sensitivity 
is estimated by adding gravitational wave signals into ei- 



'h+(t) 
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cos(2nf t) 


exp 




h x (t)_ 


~d 


_sin(27r/ i)_ 
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ther real or simulated detector noise; approximately 10 4 
such "injections" are performed spread over 2 decades in 
distance. We use as our model signal a circularly polar- 
ized Gaussian-modulated sinusoid with central frequency 
fo = 150 Hz and Q = 9, 



(20) 



Here A is an arbitrary scaling factor and d the distance 
to the source. This is a standard waveform for evalua ting 
the sensitivity of triggered GWB searches @, H [H, 0, 

For an ideal matched filter search in Gaussian noise 
d 50 % should correspond to the distance at which the me- 
dian of p c for that signal model crosses a certain thresh- 
old. This threshold depends on the number of degrees of 
freedom of the filtering template and the effective number 
of independent times. For the search considered here, the 
number of independent trials is of the order of 3 x 10 5 : 
the total time-frequency volume in the on-source win- 
dow. The false alarm probability per trial is thus ap- 
proximately 3 x 10~ 8 . The number of degrees of freedom 
is 2 (real and imaginary parts of the data) times the num- 
ber of time-frequency pixels in a Gaussian noise cluster 
(typically in the 4—6 range). Hence using a \ 2 distribu- 
tion for Gaussian noise we obtain an expected threshold 
on p c in the 7.1—7.7 range. 

However, for an analysis of real data, p c is not neces- 
sarily a good figure of merit, as a signal needs also to pass 
coherent consistency tests to be distinguished from non- 
Gaussian noise transients. Heuristically, a signal needs 
a sufficient E — I difference that Gaussian noise fluctu- 
ations will not destroy the signal consistency as seen by 
the coherent tests described in Sec. HVl Empirically we 
find that the penalized coherent SNR 
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(21) 
(22) 



where A^et is the number of detectors in the network, is 
a good figure of merit for signals in non-Gaussian noise 
from real gravitational wave detectors [13] ■ The normal- 
ization factor (Adet + l)/(A r d e t — 1) is chosen so that the 
penalty factor A c /p c has a maximum value of 1. Note 
that the magnitudes of the right and left-handed sensi- 
tivity vectors are equal 



l/ Q °IH/ a °l = (/+ 2 + /* 2 ) 



1/2 



(23) 



so it does not matter which is used to evaluate Eq. 

The distribution on the sky of the penalty factor for 
different networks is shown in Fig. [5J The maximum of 
the factor for any network is attained when the sensitivity 
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FIG. 3. Values of p c and A c corresponding to the detec- 
tion sensitivity dso% found when analyzing real and simulated 
Gaussian noise. Each marker represents a complete end-to- 
end analysis for one of the networks and sky positions indi- 
cated in Fig. [2] Plus marks are for real noise from the H1L1 
network, cross marks for real noise from the L1V1 network, 
star marks for real noise from the H1L1V1 network, and cir- 
cles are for simulated Gaussian noise from the H1L1V1 net- 
work. The analyzed sky position are shown as crosses on 
Fig.0 
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FIG. 4. Cumulative distribution over the sky of the sensitive 
distance dso% assuming that VI is a factor 2 less sensitive 
than HI and LI. The lines show the distribution when d 50 % is 
assumed to be given by a threshold on X c (real data analysis): 
blue dashed line for the H1L1 network, red dotted line for the 
L1V1 network and black solid line for the H1L1V1 network. 
The marks show the distribution when d 50 % is assumed to 
be given by a threshold on p c (Gaussian data analysis): blue 
pluses for the H1L1 network, red crosses for L1V1 network 
and black stars for the H1L1V1 network. 



of all detectors in the network is equal, that is when all 
\f®\ are equal. For the 2 detector networks the zero 
points of the penalty factor correspond to the blind spots 
of one of the two detectors in the network. 

The test sky positions used in the comparison are also 
shown in Fig. [2] They are chosen along a line which 
samples a wide range of penalty factors, in order to dis- 
tinguish A c and p c . We use a larger sample of test sky 
positions for the LI VI network, as this is the network for 
which the penalty factor effect is most important. 

For each of the trigger sky positions and networks in- 
dicated in Fig. [5] we perform a complete analysis using 
both simulated Gaussian noise and real data. In each 
case we compute the distance sensitivity d$Q% at which 
the source model (|20| is detectable with 50% probability 
at a fixed false alarm probability of 1%. The values of 
p c and A c corresponding to this d 50 % are then computed 
for each test sky position; see Fig. |3] For simulated noise 
we obtain the expected result that corresponds to a 
threshold on p c , here equal to ~ 7.4, which falls into the 
expected range of 7.1 — 7.7. However, for real noise the 
obtained value of p c is spread over a much larger range of 
7.5 — 30, whereas A c fluctuates by only 5% around ~ 7.5. 
Two conclusions can be drawn from these results: 

1. A c is a good predictor of the analysis sensitivity in 
real noise as a function of sky position for a given 
signal model; 

2. A sensitivity as good as in Gaussian noise (A c ~ p c ) 



can be attained for sky positions where the penalty 
factor is close to 1. This occurs where all detectors 
have comparable sensitivity, which corresponds to 
the white areas on Fig. [2] 

B. Distribution of sensitive distance 

The good performance of the penalized coherent SNR 
A c in predicting the sensitivity of real data analysis allows 
us to study analytically the sensitivity sky dependence, 
and to compare it with the Gaussian noise case which is 
given by p c . 

Fig. 2] shows the cumulative distribution of the sensi- 
tive distance £^50% assuming that VI is a factor 2 less 
sensitive than HI and LI. We consider two cases: a 
Gaussian noise analysis where the sensitivity is given by 
a threshold on p c and a real noise analysis where the 
sensitivity is predicted by the same threshold but on A e . 

For the H1L1 network the real and Gaussian noise sen- 
sitivities are very close; our analysis of real data is only 
a few percent less sensitive than the ideal Gaussian case. 
However for networks including Virgo, especially for the 
L1V1 network of two non-aligned detectors, the sensitiv- 
ity is as much as 20% lower with real data than with 
ideal Gaussian noise. This is an expected effect of each 
detector having maximum antenna response near the null 
response of the other detector, which limits the spurious 
transient noise rejection methods as they rely on the sig- 
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FIG. 5. Fractional value of d 5 o% for real detector noise as 
a function of the number of analyzed sky positions relative 
to the expectation based on the typical threshold A c = 7.5, 
which we obtained in Fig. O Blue plus marks denote results 
for the H1L1 network, red cross marks for the L1V1 network 
and black circle marks for the H1L1V1 network. The 6 marks 
around 100 sky positions correspond to search sensitivity over 
the sky patches represented as circles on Fig. [2] The 2 marks 
with sky position number > 1000 correspond to search sensi- 
tivity over the full sky. 



nal being visible above the Gaussian part of the noise in 
at least two detectors. 

Interestingly for 30% of the sky the sensitivity of the 
H1L1V1 network is slightly worse than that of the H1L1 
network for the real data case. The sensitivity loss oc- 
curs for areas of the sky where the antenna patterns are 
optimal for the HI and LI detectors, and relatively poor 
for VI. Hence for these sky regions Virgo brings addi- 
tional non-stationary noise but only a very small increase 
in the total gravitational wave signal. We note that the 
sensitivity difference is less than 10%, and is a conse- 
quence of the non-optimality of our analysis of real data 
(the optimal procedure for real data is not known). To 
verify this prediction, we repeat the full end-to-end anal- 
ysis of a dozen sky positions in those regions for both 
the H1L1 and H1L1V1 network. The comparison of the 
obtained d 50 % confirms that the H1L1V1 network is less 
sensitive in those sky regions by up to 10%. This indi- 
cates that care should be taken when selecting detectors 
for the analysis of real data in mis-aligned networks. 



C. Large sky area analysis 

In principle, searching over a large sky area should 
lower the sensitivity, due to the trials factor incurred 
from repeating the analysis over a grid of sky positions. 
We have seen that the penalized coherent SNR A c is a 
good figure of merit for the sensitivity of X-Pipeline as 



a function of a single sky position. We therefore expect 
that for large sky areas the 50% efficiency distance d 50 % 
should correspond to a threshold on the median value of 
A c over the analyzed sky area, where the median takes 
into account the prior pem(^) on the true source sky 
position. Due to the trials factor, the threshold on this 
median A c will typically be slightly higher than in the 
single sky position case, and the ratio between the two 
yields an estimate of the sensitivity loss due to a large 
sky area analysis. 

To assess this sensitivity loss we repeat the analysis us- 
ing typical sky location uncertainties for the GBM instru- 
ment on Fermi: a 5° statistical error and two-component 
systematic error as described in Ref. This results in 
a search grid of ^700 square degrees in area. Due to lim- 
ited computational resources only a small number of sky 
regions and network combinations are used; the analysed 
sky regions are shown as circles on Fig. [2] For complete- 
ness, we also perform a full-sky analysis for the H1L1 and 
L1V1 networks. 

The resulting sensitive distances relative to the 

expected value for a single sky point analysis are shown 
in Fig. El We find that the performance loss of GBM- 
type error regions compared to precisely-localized exter- 
nal triggers is less than 10%. A complete lack of sky po- 
sition information (requiring an all-sky search) decreases 
the sensitivity by ~ 20%. Hence the availability of exter- 
nal triggers which are well localized on the sky (< 1 deg 2 ) 
improves the sensitivity by up to 20%. This is in addition 
to the sensitivity improvement resulting from the known 
time of the trigger, which reduces the trials factor from 
the length of data to be searched [38j . 



VI. CONCLUSION 

We have studied how the sensitivity of a search for 
GWBs performed by X-Pipeline depends on the sky re- 
gion specified by an external astrophysical trigger. Two 
aspects of the sky region affect the search sensitivity: the 
magnitudes of the antenna patterns of the various grav- 
itational wave detectors in the network over this region; 
and the area of the region, which affects the size of the 
parameter space of the search. 

For the case of Gaussian background noise we have ob- 
tained the result expected for an optimal analysis, namely 
the sensitivity is given by a threshold on the coherent 
SNR p c , and this threshold falls into the range predicted 
by a X 2 distribution given the number of independent 
trials in the search. 

For real data we introduce a penalized coherent SNR 
A c , which proves to be a good predictor of the search 
sensitivity in real non-Gaussian noise. It is expressed 
as the coherent SNR times a penalty factor; this penalty 
factor is equal to 1 (no penalty) if all detectors have equal 
sensitivity for a given sky position, and to if only one 
of the detectors is sensitive. For regions of the sky where 
the penalty factor is equal to 1, the sensitivity in real 



9 



noise is as good as in the Gaussian noise case. 



all-time searches for gravitational waves. 



The penalized coherent SNR A c allows us to separate 
the effect of antenna patterns changing over the sky from 
the effect of a search parameter space increase due to a 
large search sky region. We find that trigger informa- 
tion that allows us to restrict the search to a single sky 
position increases the sensitivity by ~ 20% compared to 
searching over the whole sky, and by ~ 10% compared to 
searching over error boxes of a few hundred square de- 
grees. This addresses part of the long-standing question 
on how externally triggered searches relate to all-sky and 
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